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Abstract 

Using the updated proton and helium fluxes just released by the Ams- 
02 experiment we reevaluate the secondary astrophysical antiproton 
to proton ratio and its uncertainties, and compare it with the ratio 
preliminarly reported by Ams-02. We hud no unambiguous evidence 
for a signihcant excess with respect to expectations. Yet, some prefer¬ 
ence for a flatter energy dependence of the diffusion coefficient (with 
respect to the Med benchmark often used in the literature) starts to 
emerge. Also, we provide a hrst assessment of the room left for ex¬ 
otic components such as Galactic Dark Matter annihilation or decay, 
deriving new stringent constraints. 
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1 Introduction 

Since decades, the antiproton [p) component in cosmic rays (CR’s) has been recognized as 
an important messenger for energetic phenomena of astrophysical, cosmological and particle 
physics nature (see for instance [1, 2, 3]). In modern times, antiprotons have often been argued 
to be an important diagnostic tool for CR sources and propagation properties, and constitute 
one of the prime channels for indirect searches of Dark Matter (DM) [4, 5], which so far has 
only been detected gravitationally. In DM annihilation (or decay) modes, antiprotons can result 
either from the hadronization of the primary quarks or gauge bosons or through electroweak 
radiation for leptonic channels. 

The Alpha Magnetic Spectrometer (Ams-02) onboard the International Space Station (ISS), 
is the most advanced detector for such indirect DM searches via charged CR flux measurements. 
The positron fraction has been published earlier [6, 7], confirming the rise at energies above 
10 GeV detected previously by Pamela [8, 9] and Fermi [10]. The sum of electrons and 
positrons [11] as well as their separate fluxes [12] have also been published, thus drawing a 
coherent and extremely precise picture of the lepton components of CR’s. Despite the fact 
that DM interpretations of the positron and, more generally, leptonic ‘excesses’ have been 
attempted (for a review see [13]), even before the advent of Ams-02 it had been recognized 
that explanations involving astrophysical sources were both viable and favoured (for a review 
see [14]), a conclusion reinforced by updated analyses (see [15, 16], and references therein, for 
recent assessments). 

In this paper, we will instead focus on CR antiprotons. Until now, the so-called secondary 
antiprotons (originating from collisions of CR primaries with the interstellar material) have been 
shown to account for the bulk of the measured flux [18], thus allowing to derive constraints on 
the DM parameter space and to compute expected sensitivities, respectively based on updated 
Pamela data [19] and projected Ams-02 data (see e.g. [20, 21, 22, 23, 24, 25, 26, 27]). The Ams- 
02 Collaboration has now presented its preliminary measurements of the p/p ratio [29], with an 
improved statistical precision and energy range extending to 450 GeV. It is therefore crucial and 
timely to re-examine the situation and update existing results. In addition, Ams-02 has pub¬ 
lished the measurement of the proton (p) spectrum [28] and presented the measurement of the 
helium (a) one [29], in qualitative agreement with the previous determinations by Pamela [30], 
but now with unprecedented precision and detail. This is important for our purposes since the p 
and a spectra are crucial input ingredients in the computation of the secondary antiproton flux, 
which is the minimal astrophysical antiproton background for indirect DM searches, as we will 
remind later. Hence, with the release of these exquisitely precise datasets, Ams-02 provides a 
coherent, high-statistics—albeit preliminary—picture in the hadronic component of CR’s too, 
allowing for a scrutiny of possible exotic contributions. 
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However, the reach of any search for exotic physics is limited by the astrophysical uncer¬ 
tainties affecting the production and the propagation processes of cosmic antiprotons in the 
Galaxy and in the solar system. Indeed, while the basic processes involved in the production 
and propagation of CR antiprotons are rather well understood, the detailed parameters entering 
in such processes are far from being well determined. The p production, propagation and Solar 
modulation uncertainties can have a large impact on both the astrophysical and (in particu¬ 
lar) the DM signal. Some sensible ranges for these parameters can and must be determined 
by studying ordinary CR fluxes like the ratio of Boron to Carbon (B/C ratio), which surely 
have a non-exotic origin. Indeed, in this way the traditional Min-Med-Max schemes [31] are 
determined, and plausible ranges for the force field parameter of solar modulation (the so-called 
Fisk potential) are identified. However, these ranges are based on past CR data and are not 
necessarily guaranteed to work in describing the current status. We anticipate that this is what 
we will hnd in some cases discussed below. For instance, a string of recent papers, based on 
synchrotron radio emission [32, 33, 34, 35] but also on positrons [15, 36] and somewhat also 
on gamma rays [37], hnds that the thin halo predicted by Min is seriously disfavored. More 
generally, looking for DM on top of inadequate schemes can lead to non-robust or even wrong 
conclusions. Hence, one of the most crucial issues in the field is to update the uncertainty ranges 
of ordinary astrophysics in view of the more recent and precise experimental results, in order 
to build the DM search on a more solid basis. This will be possible after a careful analysis of 
accurate secondary over primary data like the B/C ratio, soon to be published by Ams-02 (and 
possibly other experiments), and provided that theoretical uncertainties will be under better 
control [38] . For the time being, the search for DM signatures has to be pursued with the utmost 
care. 

Within this broad context, the purpose of this paper is twofold: 1) based on existing prop¬ 
agation models, derive the state-of-the-art astrophysical antiproton background, carefully ap¬ 
praising the related uncertainties; 2) on the basis of such background and fully taking into 
account such uncertainties, assess what can be said on the room left for a DM signal, and what 
can not. 

The rest of this paper is organized as follows. In Sec. 2, we remind how the computation of the 
astrophysical antiproton background proceeds, we detail its uncertainties and we compare the 
result with the measured p/p. In Sec. 3 we introduce the DM contribution to p/p and we derive 
constraints on the DM annihilation cross section or decay rate, for several annihilation/decay 
channels and under different DM and astrophysical configurations. Finally in Sec. 4, we conclude 
with a few final comments. 


2 Re-evaluation of the astrophysical antiproton background 

The secondary astrophysical antiproton background ^ is produced in collisions of the CR high 
energy protons and helium nuclei on the interstellar medium, mainly constituted of hydrogen 
and helium, the contributions of heavier nuclei in both projectiles and targets being a few 
percent correction. The locally measured flux is the result of the diffuse production in the 
Galactic environment and the subsequent propagation of the antiprotons to the location of the 
Earth. Hence, the main ingredients of the computation for the ‘(secondary) astrophysical p 

^In some models one can also have a primary astrophysical source of background antiprotons, i.e. a significant 
antiproton population participating to the acceleration process, see e.g. [39]. 
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source term’ are: i) the injection p and a primary fluxes from Galactic sources, ii) the collision 
cross sections, iii) the propagation details. While we refer to [21, 40] and reference therein for 
a detailed discussion of all the aspects of the computation, here we just highlight the points of 
novelty. 

For the p and a spectra needed in i), as mentioned above we use the data that have just been 
released by Ams-02 [28, 29]. The spectra are measured up to a rigidity of 1.8 and 3 TV for p and 
a nuclei, respectively, and, as already reported by the Pamela Collaboration [30] , they cannot 
be described by a single power law: a hardening at energies higher than ~300 GV is observed 
for both. At the practical level, we perform our own fits of the Ams-02 data points. The value 
of the Fisk potential which gives the best our fits is (f)p = 0.62 GV, the upper bound of 

the interval sets in [28]. The values of the best-fit parameters are reported in appendix. The 
uncertainties on the slope of the p and a spectra at high energies, induce an uncertainty 

band on the predicted astrophysical p/p ratio. In fig. 1, top left panel we show the result of 
our computation of the ratio with such uncertainty band. For the distribution of the sources of 
primary CR p and a, which can be determined from pulsar and supernova remnant surveys, we 
use the parameterization of [42], slightly modified as in [43]. 

For the production processes we need the cross sections UpH-s-pX, UpHe-s-pX, UcH-s-px, UaHe^-pX, 
where the first index refers to the impingent primary CR while the second one to the target 
interstellar material. For UpH we use the new parameterization recently proposed by [44] , instead 
of the traditional fitting relations given in [48, 49]. For the cross sections of the other reactions 
we use the prescription of [40], to which we refer the interested reader. We just remind that 
for the cross section values that we adopt the pH reaction dominates, providing 60% to 65% 
of the total p flux depending on the energy, while pHe and aH reactions yield 32 to 37%, and 
the reaction aHe contributes less than 3%. Another element which has only recently been 
appreciated is related to the contribution of antineutron production: on the basis of isospin 
symmetry, one would consider the production cross section for antineutrons (e.g. CTpH-s-nx and 
the others) as equal to those for antiprotons; the antineutrons then rapidly decay and provide an 
exact factor of 2 in the p flux. However, as pointed out in [44, 50] and as already implemented 
in [21], it may be that this naive scaling does not apply and that the antineutron cross section is 
larger by up to 50% with respect to the p one. Assessing uncertainties for reactions involving He 
is even more challenging, since no data are present, and predictions are based on semi-empirical 
nuclear models calibrated on data involving either protons or heavier nuclei (see [51]). For 
sure, uncertainties involving these reactions are at least as large in percentage as the one of the 
pH reaction, an assumption we will do in the following. More conservative assumptions would 
only make the error larger, and strengthen our main conclusion on the level of agreement of 
the data with a purely secondary antiproton flux. All these cumulated effects contribute to an 
uncertainty band for the astrophysical p/p ratio which is represented in fig. 8 of [44] and which 
we will adopt: it varies from about 20% to at most 50% (at large energies and in the most 
conservative conditions). In fig. 1, top right panel, we show our prediction for the p/p ratio 
with this uncertainty envelope. 

Once produced, antiprotons have to propagate in the local Galactic environment before they 
are collected at Earth. We deal with this process in the usual way, by solving semi-analytically 

^Each of the panels of the figure has to assume a choice for the uncertainties presented in the other panels. 
E.g. the first panel assumes definite values for the collision cross sections, a model for p propagation and a value 
for the Fisk potential. They are always chosen to be the central values, e.g Med, the fiducial cross section and 
0.62 GV for this example. 
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Figure 1: Illustration of the individual partial uncertainties for secondary antiprotons. The 
colored bands represent the uncertainties on the input p and a fluxes (upper left panel), p production 
cross sections in the interstellar medium (upper right panel), Galactic propagation (lower left panel) 
and Solar modulation (lower right panel). 


the full transport equation for a charged species in a 2D cylindrical ‘thick halo’ model of the 
Galaxy. We do not reproduce the full treatment here (we refer again to [21] for a self-contained 
description and to [45, 46, 47, 55, 31, 40] for all the relevant details) but point out that we do 
include all the relevant processes. In particular, we take into account p annihilation, energy 
losses, ‘tertiary production’, and diffusive reacceleration. Besides these effects, the propagation 
parameters governing diffusion and convection are as usual codified in the Min, Med and Max 
sets [31], which are by dehnition those that minimize or maximize a hypothetical primary, DM p 
flux at Earth. Note that these have not (yet) been revised on the light of recent secondary data 
like the preliminary B/C ratio of Ams-02, as discussed in the introduction, so the viability of 
these predictions for the p/p ratio (which extends for instance to higher energies) is not trivially 
expected to hold. In Eg. 1, lower left panel, we show the impact of the propagation uncertainty. 
The curves which are labelled Min, Med and Max represent the modification which occurs by 
choosing these standard sets. The shaded yellow area envelops the results obtained by sampling 
more widely the propagation parameter space that has been shown in [55] to be compatible 
with the B/C ratio and finding the values that minimize and maximize the secondary, rather 
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Figure 2: The combined total uncertainty on the predicted secondary p/p ratio, superim¬ 
posed to the older Pamela data [53] and the new Ams-02 data. The curve labelled ‘fiducial’ assumes 
the reference values for the different contributions to the uncertainties: best fit proton and helium 
fluxes, central values for the cross sections, Med propagation and central value for the Fisk potential. 
We stress however that the whole uncertainty band can be spanned within the errors. 


than primary, p/p flux. Notice that the shaded yellow area does not coincide with the Min- 
Med-Max envelope (see in particular between 50 and 100 GeV): this is not surprising, as it 
just reflects the fact that the choices of the parameters which minimize and maximize the p/p 
secondaries are slightly different from those of the primaries. However, the discrepancy is not 
very large. We also notice for completeness that an additional source of uncertainty affects the 
energy loss processes. Among these, the most relevant ones are the energy distribution in the 
outcome of inelastic but non-annihilating interactions or elastic scatterings to the extent they 
do not fully peak in the forward direction, as commonly assumed [55]. Although no detailed 
assessment of these uncertainties exists in the literature, they should affect only the sub-GeV 
energy range, where however experimental errors are significantly larger, and which lies outside 
the main domain of interest of this article. 

Finally, p’s have to penetrate into the heliosphere, where they are subject to the phenomenon 
of Solar modulation (abbreviated with ‘SMod’ when needed in the following figures). We de¬ 
scribe this process in the usual force held approximation [52], parameterized by the Fisk po¬ 
tential (pF, expressed in GV. As already mentioned in the introduction, the value taken by fip 
is uncertain, as it depends on several complex parameters of the Solar activity and therefore 
ultimately on the epoch of observation. In order to be conservative, we let fip vary in a wide 
interval roughly centered around the value of the fixed Fisk potential for protons ffp (analo¬ 
gously to what done in [25], approach ‘B’). Namely, fip = [0.3,1.0] GV ~ 0^ ± 50% 0^. In 
fig. 1, bottom right panel, we show the computation of the ratio with the uncertainties related 
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to the values of the Fisk potential in the considered interval. Notice hnally that the force held 
approximation, even if ‘improved’ by our allowing for different Fisk potentials for protons and 
antiprotons, remains indeed an “effective” description of a complicated phenomenon. Possible 
departures from it could introduce further uncertainties on the predicted p/p, which we are not 
including. However it has been shown in the past [24] that the approximation grasps quite well 
the main features of the process, so that we are conhdent that our procedure is conservative 
enough. 

Fig. 2 constitutes our summary and best determination of the astrophysical p/p ratio and 
its combined uncertainties, compared to the new (preliminary) Ams-02 data. The crucial 
observation is that the astrophysical hux, with its cumulated uncertainties, can reasonably well 
explain the new datapoints. Thus, our hrst —and arguably most important— conclusion is 
that, contrarily to the leptonic case, there is no clear antiproton excess that can be identified 
in the first place, and thus, at this stage, no real need for primary sources. This also means 
that, at least qualitatively, one expects a limited room left for exotic components, such as DM. 
Indeed in the following section we will proceed to compute the constraints on it. 

However, before we can do so, we have to identify specihc sets of astrophysical parameters to 
describe the background, as discussed in the introduction. We £x in turn Min, Med and Max 
and we vary the Solar modulation potential in the given interval. We model the uncertainties 
of the production cross sections term by allowing a renormalization of the background with an 
energy dependence and an amplitude A as dictated by the analysis presented above (namely, 
an uncertainty modulated as the pink band of fig. 1). With this strategy, we look for the best 
fitting values of the amplitude A and of the potential fip and we trace the corresponding p/p 
spectra. In concrete terms, for each propagation model, we minimize the chi-square Xq(H, 
with respect to the Ams-02 data and hence determine the best £t amplitude Aq and Fisk 
potential We show in £g. 3 the different cases. 

Even within the limitations of the data like those we are dealing with (namely their prelimi¬ 
nary nature, their errors only partially accounted for and the partial collection time with respect 
to the full lifetime of the experiment), we can see that the Min propagation scheme predicts 
an astrophysical background that can not reproduce the new p/p data points above 30 GeV. 
The Med scheme provides a barely decent fit (still good up to ~ 30 GeV but rapidly degrading 
after) while choosing Max the data can be well explained across the whole range of energies. 
We have explicitly computed the corresponding yfi to support the above statements, with the 
Min, Med, Max cases yielding 106, 58 and 41, respectively (for 28 degrees of freedom). Given 
the preliminary nature of the data, of course they have only an indicative significance. This is 
our second conclusion: the preliminary p / p Ams-02 data seem to prefer a model, such as Max, 
characterized by a relatively mild energy dependence of the diffusion coefficient at high energies. 
Although it is too early to draw strong conclusions, this is an interesting observation and it goes 
in the same direction as the preference displayed by the preliminary B/G Ams-02 data [38] 

It would of course be tempting to interpret the room left in the Min and Med cases at 
large energies as an exotic contribution from DM. However we insist that this would be a wrong 
deduction in two respects: as long as a model within the uncertainties can £t the data, failure of 
other models just means a better selection of the background rather than evidence for an extra 

At is also backed by the results recently reported in [56]— appeared after the first version of this paper in 
pre-print form—based on fits to Pamela B/C, p and He data. 
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component; in any case, a new assessment of the viable propagation parameter space would be 
needed before any conclusion is drawn. 


3 Updated constraints on Dark Matter 

Primary antiprotons could originate from DM annihilations, or decays, in each point of the 
Galactic halo. They then propagate to the Earth subject to the same mechanisms discussed in 
the previous section, which are in particular described by the canonical sets of parameters Min- 
Med-Max. Concretely, we obtain the p fluxes at Earth (post-propagation) from the numerical 
products provided in [54], version 4. Notice that these include the subtle effects of energy 
losses, tertiaries and diffusive reacceleration which, as discussed at length in [21], are important 
to reach a detailed prediction. 
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We consider four primary annihilation (or decay) channels: DM DM —)■ bb,W~^W~, iJ.~ 
and 77 . These, for all practical purposes, cover very well the range of possible spectra. Indeed, 
annihilation (or decay) into ti or hh (with h the Higgs boson) would give spectra practically 
indistinguishable from those from DM DM —)■ bb, and ZZ from those of W~^W~. The 
channel represents leptonic channels, in which a small yield of antiprotons is obtained thanks 
to electroweak corrections (namely, the radiation from the hnal state leptons of a weak boson 
which decays hadronically). Similarly, the 77 channel produces some subdominant p flux via 
electromagnetic corrections 

We also consider two representative DM Galactic prohles: Einasto and Burkert, with the 
precise functional forms and dehnitions of the parameters as in [54], The former possesses a 
peaked distribution towards the Galactic center and hence typically results in a more abundant 
yield of antiprotons with respect to the latter, which features a core in the inner few kpc. 

We remind that, in section 2, we have obtained the re-evaluated astrophysical background 
fluxes and their uncertainties. In particular, we have computed the fluxes for the Min, Med 
and Max cases, displayed in Eg. 3. Armed with those and with the fluxes from DM as just 
presented, we can now compute the constraints in the usual planes ‘mass uidm vs. thermally 
averaged annihilation cross section (cm)’ or ‘mass ttidm vs. decay rate T’. We refer to [21] for 
a detailed discussion of the practical procedure, of which we just repeat here the main lines. 
Having hxed a propagation model, for a given DM mass tudm and annihilation cross-section 
{av) (or decay rate T) we add the DM signal to the secondary background. The total flux is 
then 

*^>tot(nrDM, (cm), A, (jxp) = <hbkg(^, 0 f) + *^>DM(inDM, (cm), ((>f) (1) 

and we again End the best Et amplitude and Fisk potential. Finally, we solve the following 
equation in (cm) 

(cm). A, (j)F) - Xo = 4, (2) 

(where Xo is the minimum chi-squared of the background-only case as computed in the previous 
section) in order to obtain the exclusion contour. We reproduce this for each mass point. 

On the basis of the discussion in the previous section, it makes sense to derive constraints 
only within the propagation schemes that provide a decent explanation of the background. Max 
is the favored scheme. Med provides overall a worse but still reasonable Et to the data, so that 
we will employ it. In addition, (see Eg. 3, middle panel) at small energies (T < 30 GeV) its Et 
is good, thus meaningful constraints on relatively light DM (tudm ^ 300 GeV) can be derived. 
We discard instead the Min case. 

The results that we obtain with this strategy are presented in Eg. 4 for the DM annihilation 
case and in Eg. 5 for the DM decay case. In the left panels we Ex a benchmark DM proEle 
(Einasto) and the Med propagation model, and show the constraints for the diEerent particle 
physics channels introduced above. We see for example that the thermal annihilation cross 
section (cm) = 3 ■ 10“^® cm^/s is now touched by the exclusion line for tudm ~ 150 GeV for the 
bb channel. In the right panels we explore the impact of changing the propagation parameters or 
the DM distribution. As already highlighted several times in the literature, the efiect is sizable 
and can reach a factor of up to an order of magnitude. For instance, the previously quoted 
limit for the mass of a thermal relic can vary between 90 and 250 GeV for the range of models 
explored here. Of course, as Max maximizes by deEnition the DM p yield, its constraints are 

“^For simplicity, we consider only the production of p from the final state. In principle, in this channel, 
additional hadronic production is possible from the states mediating the process of DM annihilation into photons. 
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Annihilation constraints from p j p Astrophysical uncertainties on the constraints 




Figure 4: Annihilating DM: current constraints. Left Panel: current constraints from the 
antiproton to proton ratio measurements by Ams-02, for different annihilation channels. The areas 
above the curves are excluded. Right Panel: illustration of the impact of DM-related astrophysical 
uncertainties: the constraint for the bh channel spans the shaded band when varying the propagation 
parameters (dashed lines) or the halo profiles (solid lines). Notice that in the Min case the analysis is 
not sensible, hence not shown here (see text for details). 


Decay constraints from p / p Astrophysical uncertainties on the constraints 




Figure 5: Decaying DM: current constraints. Left Panel: current constraints from the antiproton 
to proton ratio measurements by Ams-02, for different decay channels. The areas below the curves 
are excluded. Right Panel: illustration of the impact of DM-related astrophysical uncertainties: the 
constraint for the bh channel spans the shaded band when varying the propagation parameters (dashed 
lines) or the halo profiles (solid lines). Notice that in the Min case the analysis is not sensible, hence 
not shown here (see text for details). 
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much stronger than those of the Med case. Turning the argument around, if the preference for 
MAX-like propagation schemes hinted at by preliminary Ams- 02 data is confirmed, Ams- 02 
itself has the unprecedented possibility to exclude tudm ^ 250 GeV for thermal annihilation 
cross section in the bb channel. 


4 Final remarks 

In the light of the new p flux published by Ams-02 and the preliminary Ams-02 results pre¬ 
sented on the a flux as well as the p/p ratio, and using the new results of the p production cross 
sections, we have re-evaluated the secondary astrophysical predictions for the p/p ratio. We 
have accounted for the different sources of uncertainties: namely on the injection fluxes, on the 
production cross sections, on the propagation process and those connected to Solar modulation. 
Our first and main result is that there is no unambiguous antiproton excess that can he identified 
in the first place, and thus, at this stage, no real need for primary sources of antiprotons. Within 
errors, secondary astrophysical production alone can account for the data. This conclusion is 
highly non-trivial, since we relied on updates of existing propagation schemes, which were not 
necessarily expected to work in the high precision and extended energy regime made accessible 
by Ams-02. Adopting a more conservative treatment of the uncertainties of antiproton produc¬ 
tion cross-sections involving He as either target or projectile nuclei would clearly reinforce this 
conclusion. 

Next, we enter in the merit of which propagation schemes do account for the data, taking 
into account the other uncertainties. We find that the data seem to prefer a model, such as 
Max, characterized by a relatively mild energy dependence of the diffusion coefficient at high 
energies. If confirmed, this would go in the same direction as other indications already obtained 
in different channels, as discussed above. 

Finally, an important application concerns updated constraints on DM: within the framework 
of the propagation schemes that it is sensible to use, we derive bounds that are more stringent 
by about one order of magnitude with respect to the previous ones [20, 21] (based on Pamela 
data). 

Of course, this analysis is very preliminary and there is still room for improvements. First 
and foremost, the release of the final p/p measurement with systematic and statistical errors 
fully accounted for. Yet, even a preliminary analysis allows to show that antiprotons confirm 
themselves as a very powerful probe for CR physics and for DM in particular. Actually, consid¬ 
ering the puzzling excesses (with respect to the originally predicted astrophysical background) 
of undetermined origin in the electron and positron fluxes, considering the complicated back¬ 
ground of most gamma-ray searches and considering the challenges of neutrino detection, p’s 
might arguably still be the most promising avenue in DM indirect searches, since improving 
the knowledge of the background is relatively easier than for other channels and so perhaps 
seeing the emergence of a clear signal is possible. In this respect, the Ams-02 experiment can 
play a crucial role. So far it has essentially confirmed the results of previous experiments (most 
notably Pamela), but it has done so with an impressively improved accuracy: the qualitative 
picture in DM indirect searches has been left largely unchanged by it, but Ams-02 has allowed 
improved pinning down of the parameters and tightening of the constraints. In this context, 
while follow-up releases of antiproton data (e.g. pure fluxes, extended energy range or enlarged 
statistics) will obviously be welcome, it is urgent to address first one of the main current limi¬ 
tations in the field of charged CRs, namely the determination of the propagation parameters. 
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In this respect, analyzing the upcoming reliable and accurate light nuclei measurements from 
Ams-02 will provide the community with a very powerful leverage for any search of exotics in 
CR’s. At that point it will be possible to assess whether or not excesses are present in antiproton 
data, for instance if the current small deficit increases in significance (although identifying their 
origin will remain very challenging [57]). 


Appendix: Parametrization of primary proton and helium 
fluxes 


To fit the Ams-02 p and a fluxes we used the following rigidity dependent function (in particles 


m-^s-^sr-^GV-^l 


$ = G • (1 - 


R 


R -|- (j)p 


{R + M 


1 + 


R (j)p 
Rb 


A7 


(3) 


We proceed in two steps. First 7 , Ay, Rp, s are fixed using the high energy part {R > 45 GV) 
of the spectrum. Then G, a and f3 are determined over the all energy range. The value of the 
Fisk potential which gives the best for our fits is = 0.62 GV, the upper bound of the 
interval sets in [28]. The values of the best-fit parameters are reported in table 1. 



V 

a 

G 

23566± 30 

4075± 2 

A 

-0.519 ± 0.007 

-0.163 ± 0.004 

/S 

1.21 ± 0.02 

0.41 ± 0.02 

7 

-2.849 ± 0.002 

-2.795 ± 0.009 

Rb 

355 ± 33 

284 ± 38 

Ay 

0.146 ± 0.02 

0.162 ± 0.009 

s 

0.0325± 0.0131 

0.078± 0.035 

"X^ndof 

29.02/(73-7) 

2.62/(54-7) 


Table 1; Best-fit values for p and a fluxes. 
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